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ABSTRACT 


Using the ISE£-3 ra^lo astronomy experiment data we nave identified 37 « 

interplanetary (IP) type II bursts in the period September 1978 to December 
1981. We List these events and the associated phenomena. The events are 
preceded by intense, soft X ray events with long decay times (LDEs) and type 
II and/or type IV bursts at meter wavelengths. The meter wavelength type II 
bursts are usually intense and exhibit herringbone structure. The extension of 
the herringbone structure into the Icilometer wavelength range results in the 
occurrence of a shoclc accelerated (SA) event. The SA event is an important 
diagnostic for the presence of a strong shock and particle acceleration. The 
majority of the interplanetary type II bursts are associated with energetic 
particle events. Our results support other studies which indicate that 
energetic solar particles detrcced at 1 A.U. are generated by shock 
acceleration. From a prelimir<ary analysis of the available data there appears 
to be a high correlation with white light coronal transients. The transients 
are fast i.e. velocities greater than 500 km/sec. 
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I. 

ABSTRACT 

Ualng the 12££-3 radio astronomy Experiment data we have identified 37 
interplanetary (I?) type II bursts in the period September 1978 to December 
1981. Ue list these events and the associated phenomena. The events are 
preceded by intense, soft X ray events with long decay times (LDEs) and type 
II and/or type IV bursts at meter wavelengths. ‘Ihe meter wavelength type II 
bursts are usually intense and exhibit herringbone structure. The extension of 
the herringbone structure into the Idlometer wavelength range results in the 
occurrence of a shock accelerated (SA) event. The SA event is an important 
diagnostic for the presence of a strong shock and part.'cle acceleration. The 
majority of the interplanetary type II bursts are associated with energetic 
particle events. Our results support other studies which indicate that 
energecic solar particles detected at 1 A.U. are generatec by shock 
acceleration. From a preliminary analysis of the available data there appears 
to be a high correlation with white light coronal transients. The transients 
are fhst i.e. velocities greater than 500 km/ sec. 
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1. INTilOOaCTION 

Th« typa II burst results from plasma emission generated by a shocx as it 
propagltes out through the solar corona. The frequency drift rate of the burst 
1s related to the shock's velocity. Ground-based observations of type II 
bursts pertain to coronal heights leas than about 5 solar radii. A shock 
typically takas 2 days tc r^ach the earth at 215 solar radii vhere it causes a 
sudden commencement geomagnetic storm . In situ observations of shocks provide 
information about shock properties including the velocity, but such 
observations have been made primarily by earth orbiting satellites. Thus shock 
properties in the region between 5 and 215 solar radii are not well studied. 
Satellite experiments operating at low radio frequencies can access this 
region by remotely observing Interplanetary type II bursts. The value of such 
observations are limited in part by our current understanding of the radio 
emission process. 

Two events were reported f1*om observations from IMF-6 during the previous 
period of solar maximum (Malitson et al., 1973*1976). Boischot et al. (I960) 
reported the detection of a number of events fl”om the Voyager spacecraft and 
there have been accounts of detections from the Prognoz-6 satellite (Pinter et 
al., 1962). However the IS££-3 radio astronomy experiment (Snoll et al., 

1978), which operates over a frequency range of 2 MHz to 30 kHz, has obtained 
the moat complete and detailed set of observations of interplanetary (I?) t:rpe 
II bursts thus far available. The experiment is more sensitive than previous 
experiments and because of its orbit is much less troubled by terrestrial 
kUometric radiation. In addition the experiment can observe the sun 
continuously. An Initial paper reporting the detection of 12 events was 
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prepared in 1980 and published in 1982 (Cane et al. 
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Cane et al. (1981) reported on a new class of kHometer wavelength 
bursts, the shock accelerated (SA) events, alsb related to energetic shocks. 
The SA event is the low frequency continuation of the L'erringbone structure 
associated with meter wavelength type II bursts. At kilometer wavele;:gths 
these events precede the IP typu II burst. The SA event allows the 
determination of the start of the type II event low in the corona. The 
sequence of events is illustrated in figure 1 (from Cane et al., 198I). There 
is no frequency coverage between about 20 and 2 MHz. The tine difference 
between the end of the meter wavelength type II burst and the observation of a 
type II burst at 2 MHz is of the order of 1/2 hour. However the rapid drift of 
the SA event means that the tine difference between the 2 MHz SA event and the 
start of the meter wavelength type II burst is of the order of a minute. Thus 
associations between meter wavelength and kilometer wavelength type II bursts 
can be made unambiguously. 


In this paper we present information on the IP type II events and 
associated phenomena. Since the writing of the previous paper our 
understanding of the data has greatly Improved and the sample of events has 
Increased three- fold. It is therefore timely tc provide an up-date and more 
comprehensive description. 


II. DATA ANALYSIS 

The ISEE-3 radio astronomy data shows mmerous slow drift features in the 
dynamic spectra (plots of intensity as a function of frequency and time). The 
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najority of thta« are ahort-llved with very alow drift rataa and wa have no 
axplalnatlon for their origin, k smaller number last for many hc'ira but 
because they commence in the middle of our frequency range we have no way of 
determining a iTlkely starting time at the sun. Our list of type II bursts has 
been restricted to those events which drift through the data at a rate 
consistent with known shock velocities and which are preceeded by an SA event. 
Thus we' can identify the start of the event at the sun. Whereas other events 
may be related to solar shocks we ’include only those which are clearly the 
kilometer wavelength counterpart of the meter wavelength phenocenom. 

In table 1 we list the IP type II events. As in our previous catalog 
there are two categories*. Category 1 events have been unamoiguously Identified 
with a sudden commencement. For the moat part the type II emission is 
discemable over the frequency range at which the events are best observed 
namely 500-S0 kHz and therefore the events are observed for m^y hours. The 
reason for these upper and lower frequency bounds will be discxissed later. The 
events marked with an asterisk are not as well observed because of other 
activity occurring at the same time, which limits the detectibllity of the 
burst. 

Category 2 events are those bursts which are not followed by a sudden 
commencement or only last for a few hours. Some of the events not followed by 
a sudden commencement are associated with flares far ft*om central meridian and 
the shock probably was not extensive enough in heliographic longitude to 
intercept the earth. 

There exist a number of candidate IP type II events which are not listed. 


*5 . 
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0n« fvtnt on 0«c 5* 1981 waa txcludtd bccauat than waa no ground-baaad data 
to corroborata tha praaanca of a atrong ahoclc . Uowavar a partlcla avant waa 
datactad aa waa h auddan cofflmanceoant and it ia likely that the event 
originated behind the weat limb. Another event on May* 10, 1981 waa excluded 
bieauaa tha atart of tha avant at tha aun could not be detazinined very 
accurately. It ia probable that thia event waa aaaociatad with an eaat linb 
coronal tranaient obaarved by tha P78-1 coronograph. For a number of eventa 
there waa corroborative ground-baaed data but the low frequency data waa of 
poor quality or tha event waa obaerved only over a very amall frequency range- 

Iha table liata a number of phenomena examined in conjunction with tha 
low frequency data. In general wa have uaed data publiahed in Solar 
Gaophyaical Data (SGD) and apart f1*om Culgoora dynamic apectra, have not made 
uaa of original data aata. In the main our atudy haa been reatricted to 
phenomena which occur high in the corona. We have not uaed radio data out aide 
tha meter wavelength band. Decameter wavelength information waa not uaed 
bacauae obaervatory coverage in thia region of the apectrum ia at beat limited 
and often rendered unuaeable h*<;auae of interference. 

Tha tima givan in tha aacond column ia the atart of the meter wavelength 
type II burat or, if no type II waa reported, the atart of tha SA avent. In 
tha lattar caaa the time ia encloaed in bracketa. 

(i) Ba obaervationa 

Moat avanta hava bean aaaociatad with an Ha flare. The flare waa 
datarainad uaing tha atart of tha mater wavelength type II burat which occura 
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.within 1 ftw ainutta of tho oaxiaum in Ua (Robarta, 1959)* All tha avanta 
with good idantification wara bright flaraa of iaportanca 1 or graatar. 

Tha flara locatio'n ia alao ahown in tabla 1. Tha quaation mark danotaa an 
aaaumad bahind-tha-linb flara. Tha longitudinal diatribution of flaraa 
aaaociatad with IP typa II burata ia ahown ir tabla 2 and figura 2. 

(ii) Soft X-raya 

Tha 1-3 A aoft X-ray claaa haa baen aatimated from the daily grapha I 

I 

praaaotad in SCO. Tha majority of tha X-ray aventa are intanaa and hava dacaya 
loogar than A houra i.a thay are long duration avanta (LDEa). Tha 'T* in tha 

» 

! 

oolvmo aftar tha X ray olaaa indicataa an LDE avant waa obaarvad. Thirty- fiva | 

of tha 37 avanta wara aaan in aoft X-raya. Of the two remaining aventa, one 
occurred during an X-ray data gap and the other haa baen attributed to an 
avant behind tha waat limb. Twenty-four of the aaaociatad X-ray avanta hava \ 

long dacaya and a further 9 have dacaya between about 2 and 4 houra. Two ^ 

avanta had dacaya laaa than 2 houra. Thaaa are aaaociatad with alow, category I 

t 

2 IP typa II burata. ' 

For tha year 1931 wa have catalogued all aoft X-ray avanta whoaa 1-3 A 
claaa waa graatar than M4 and with a duration longer than 4 houra. Of Id 
avanta 10 wara aaaociatad with IP typa II burata. Six further avanta wara 
aaaociatad with SA avanta only i.a. not followed by a typa II burat. Tha 
raBaining 2 avanta oceurrad during ISEE-3 data gapa. It would appear that 
intanaa LOC X-ray avanta correlate wall with atrong ahocxa, many of which 
produce IP typa II burata. 


1 . 
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(111) M«t«r wavtltngth radio ooiiasion . 

ClgAtien of tha 37 evants ara pracadad by a oatar wavalangth type II/IV 
burst pair, 10 avants by a typa IV burst and 9 by a typa II burst. Raports of 
continuuiB ara includad undar tha classification of type IV. Tha annotations 
*W and 'P* aaan 'waak* and 'possibla* raspactivaly. Moat of tha oatar 
wavalangth typa II bursts ara clasaifiad as intanaa. tha dynaaic spactra for 
avanta also obaarvad by tha Culgoora obaarvatory axhibit cooplax behaviour 
with herringbone structure. Wa believe that all IP typa II events ara 
associated with oatar wavalangth activity. As discussed in the introduction to 
this section, the only candidate events not associated with oetar wavelength 
activity ara probably behind- tho-liob events. 

The reported occurrences of typa III bursts associated with the oeter 
wavelength type II and/or typa IV events have bean listed. The intensity class 
is given and the tine interval between the start of the oeter wavelength typa 
II burst, or tba SA event, and the reported start of the typa III activity. 

Tha annotation 'O' inplies that tha type III activity occurred during the 
oatar typa II or typa IV burst. Nine of the IP type II events are not precadad 
within 25 nins by, or associated with, ary typa III activity. Tha statistics 
do not include single bursts (i.e. type Illb) or ongoing storo activity. For 
an additional 4 avant.^ tha typa III activity connenced after the start of the 
typa II burst and nay be herringbone structure. Only one event is prceedad by 
intense type IIIG/V activity. This event occurred on July 23 19^0, during 
Culgoora observing tioe and an examination of tha data reveals the possible 
presence of two type II bursts. The first event occurred shortly after the 
type lU/V burst. The second event ccomencad a few ninutes later and is 
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th« «vtQt which oontinutd to kllootttr wavtitngtho. 

(iv) Coron&l tranalanta 

Aa ahoun in tnblt 2 thtrw w«rt l«aa IP typo II burata in 19^0 t'lan In 
1979 or 19d1. Iiona of th« 37 avanta bava baan aaaociataO with a coronal 
tranalant obaarvad by tha St9i corcnagi-apb. Convaraaly tha P-78 coronagraph; 
which bagan oparating in lata March 1979» haa obaarvad many tranaianta of 
which a numbar hava baan asaociatad with IP typa II avanta. Wa hava Indicatad 
whether a tranalant waa aaan or not with Tha 'g* indicataa a data gap. 
Thaaa gapr will be filled in aa additional data ia made available. From tha 
conpariaon to data all buc one IP typa II event haa a faat (velocity greater 
than 500 ka/aac) transient aaaoci£tad with it. For one event, oariced with a 
quaation m.vk, it ia unclear whether a transient did occur. A atudy to 
dateraine tha correlation between faat tranaienta and IP type II burata ia 
underway. 

(v) Energetic particlea 

In table 1 wa hava included tha oagnicuda of the asaociated particle 
avanta. Tha data is tha count rates frcm tha >18 Mav/n detector onboard ISE&-3 
(T.T. von Roaanvlnga, private conmunication) . Intenaity claaaea 1,2,3 
corraapond to count rataa greater than 1, 10, 100 counts/ aec respectively. 
Thirty-two of tha 37 avanta ware asaociated with energetic particle avanta 
even though oany of the flare aitea ware not wall connected. For Boat western 
avanta tha particle onset tiaa ia within an hour or two of tha start of tha 
solar activity. For acaa eastern avanta tha delay is as long aa 10 hours. 
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Thr«« «vtnta &r« oonaldtrtbly dalaytd (of thi order of 24 hours) but have be«n 
saaociated with tha solar avant bacauaa of the abaanca of other candidate 
flaraa. Thaaa 'delayed events' are marked with a gieation marie. 

In figure 3 wa show the kilooatar wavelength activity associated with 
.particle events during a period of 4 sooths. The figure suggests that larger 
particle events can be associated with I? type II bursts whereas smaller 
events can usually be associated with an SA event not followed 'by an IP type 
II burst. The association of a particle event with an SA event allows 
unambiguous identification of the associated flare, because of the positional 
information obtained with the low frequency radio experiment. 

(vi) Sudden commencements and shocks 

Sudden commencements were associated with all category 1 events by 
definition and 11 of the category 2 events. Shocks were detected at ISEE-3 
approximately 30 minutes beforri the SC. The radio astronomy experiment also 
detects the shocks in situ by the increase in the low frequency (LF) continuum 
(Hoang et. al, 1980) corresponding primarily to the increase in the ambient 
density at the spacecraft. 

(vii) Transit velocities 

The transit velocities of the shocks have been deduced from the time 
interval between the start of the event (in colmsn 2) and the arrival of the 
shock at 1 AO as determined by the sudden commencement. 

For the 22 category 1 events a mean transit tlm of approximately 2 days 
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and a a.tan transit valoolty of 840 loB/a«c is obtainad. With tht Inclusion of 
ths 11 catagory 2 avants followad by suddan coaanancamants wa find a maan 
transit valocity of 800 ka/sac. Iha velocity distribution is shovm in figura 
4. No avant had a transit valocity greater than 1100 Icn/sac (corresponding to 
a transit time lass than 1.5 Jays). 

In figura 5 via show the distribution of transit valocitlas as a function 
of flare longltxida. The dashed line shows the maan transit valocity for 6 
ranges of longitude. The distribution in longitude is reasonably uniform. TLa 
in valocity of shocks from the limbs relative to those from near centraJ. 
meridian is at most 20!i, suggesting that to a first approxlioatioc most shocks 
expand isotropically. Individual events may expand ani^iotropically such as the 
avant of Saptambar 14, 1979* This event had a transit velocity of 440 ka/sac. 
It was the last of 4 large flares from the saae active region and thus the 
shock was propagating into a very distuz*bed corona. The intensity of the event 
at kilometer wavelengths was comparable to th:.t from events with transit 
velocities near 900 ko/seC'. (We show in another paper that the intensity of an 
IP type II burst is a function of shock valocity). This suggests that radial 
above the flare site (£90) a transit valocity of the order of 900 km/ sec would 
have bean determined. 

III. 0ISCUS2I0N 

Wa have examined the correlations between I? type II bursts and other 
solar phenomena. IF type II bursts and SA avants correlate wall with energetic 
particle avants and this is consistent with theoretical (Ramaty at al., 1980) 
and observational avldanea (Gloeckler at al., 1976) which suggest that flare 
produced shocks are responsible for the .'Production of solar cosmic rays . 


ORiGii>tAL PAGE IS 
OF POOR QUALITY 


The aaaoclatlon of LDE'a with white-light tranaienta baa been ahown by 
Sheeley et al. (1975) and Kahler (1977). The aaaociation between coronal maaa 
ejection eventa and proton eventa waa ahown by Kahler et al. (1978). The 
aaaociation between particle eventa and LDE X-ray eventa la diacuaaed Jby ' 

Nonnaat et al. (1982). Since the eiirrent oodela auggeat that energetic 
partlclea are ahock accelerated all the above phenomena l.e. tranaienta, 
proton eventa and LDE X-ray eventa, ahould be aaaociated with atrong - ahocka . 
Strong ahocka are confiraed by our aaaociatlona of theae phenomena with IP 
type II burata. 

The typical atarting frequency of the fundamental of meter wavelength 
type II burata ia 7C-100 MHz (Kundu, 1965). An initial inveatigation of 
Culgoora dynamic apectra ahowa that tae atarting firequency of the meter 
wavelength btorat aaaociated with many of the IP type II eventa and the eventa 
producing an SA ovent alone, ia probably well below 70 MHz. Thia reault can be 
deduced from the obaervation of a number of SA eventa for which there waa no 
aaaociated meter wavelength type II burnt. The obaervation of an SA event at 2 
MHz impliea the preaence of a shock at coronal heigh ta below the 2 MHz plaama 
level. The type II burst ft*om which the SA event ori<;lnatea muat occur above 2 
MH?;. See Figure 1 for clarification. If no event ia detected at meter 
wavelengtha the type II burst must occur in the frequency range between 20 and 
2 MHz. Thia meana that for those SA eventa not aaaociated with a meter 
wavelength type II burst, which includes about 30^ of the eventa followed by 
an IP type II burst, the type II commenced below 20 MHz. Thia reault suggests 
that the shocks idiich survive to the IP medium are form«rfi nigh in the corona 
where presumably the absence of closed field lines facilitates their escape. 
The alternative possibility that the ahocka are formed at lower heights but 
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wr% aot producing detectable radio etaiaaion aeeaa '.eaa likely. 

The SA eventa and IP type II burata not aaaociated with a meter 
wavelength type II burat are aaaociated with L type IV burat. We auijgeat that 
the obaervation of a type IV burat ia a good indicator of the preaence of a 
ahock and that if no meter wavelength type II burat ia detected that a type II 
burat may have commenced below the loweat frequency available to ground baaed 
obaervera. 

The atandara aequence of eventa in a large flare la often illuatrated aa 
conaiating of two diatinct atagea (Wild, 2merd and Weiaa, 1963)* A group of 
intenae type III burata ia ahown to occur within a few minutea of the atart cf 
the flare. Theae are followed about five minutea later by a type II and a type 
IV burat. la can be aeen from table 1, the aequence of meter wavelength 
activity aaaociated with the IP type II eventa ia ertremely varied. Although 
the type U/IV burat pair occura about 5Qi of the time there are alao eventa 
with type II and no type IV and vice veraa. In addition the type III burata 
can coffltance before or during the type II burat and are not alwaya a aeparate 
entity. More importantly, for 9 eventa type III activity la completely 
abaent. Ihia complete abaence cf type III burata waa alao noted by Sveatka and 
Fritzova-Sveatkova (197A) who atudled the meter wavelength activity aaaociated 
with large proton eventa. The idea of a atandard aequence of eventa ia 
mialeading. 

Ue interpreted the tranait velocity diatribution aa a function of 
heliographic longitude aa indicating tont, on the average, interplanetary 
ahocka propagate iao tropically. Thia agreea with the reaulta of Chao and 



0 

ORIGINAL PAGE 13 
OF POOR QUALITY 

Lapping (1974) who found that "tha averags shock surface in the ecliptic plane 
near the earth's orbit lies on a circle centred at the sun with a radius of 1 
iU”. The fact that many of the shocks are detected at the earth and yet 
originate in regions far from central meridian, indicates the huge angular 
extent covered by such shocks. Essentially unambiguous associations can be 
made between shcclb detected at or near earth and the responsible flare region 
because of the presence of the IP type bursts and the SA events. 

Figure 2 suggests that there might he an E-W asymmetry in the location of 
flares associated with I? type II bursts. There were 15 events from sites east 
of E30 as against 9 events west of U30. We point out that because of the 
dependence of shock structure on the Archimedian spiral of the interplanetary 
magnetic field, we might expect any asymmetry to be in the eastern direction 
i.e. to favour eastern flares. The western portion of shocks are expected to 
have more highly compressed magnetic field than the eastern portion and to be 
well defined quasi-perpendicular shocks. This geometry has been invoked to 
explain' the east-west asymmetry in the magnitude of Forbush decreases 
(Bamden, 1973)* however the sample of available events is too small to 
establish the statistical significance of this result as yet. 
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IV CONCLUSION 

We have identified 37 IP type II bursts and lisled the associated 
phenomena. The following results were obtained; 

1. IP type II bursts are associated with meter wavelength type II and/or 
type IV bursts, intense LOE X-ray and energetic particle events and probably 
with coronal transients. 

2. A number of events have no associated type III activity at meter 
wavelengths . 

3. The starting 0*equencies of the associated meter wavelength type II 
bursts may be lower than average. 

The unambiguous identification of shocks detected at 1 AU with a source 
location on the sun provides the following results; 

4. The mean transit velocity of the more energetic solar shocks is 300 
km/ sec corresponding to a transit time of about 2 days. 

5. To a first approximation the shocks propagate isotropically. 


We thank T. T. von Rosenvinge for providing I2EE-3 particle data, 

R. Howard for making available P73-1 coronagraph data and R. T. Stewart for 
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Figuro 1. A schocatlc roprostntation of tht rolatlonship botwoon oetor 
wavolongth typo II activity with herringbono structure and the activity 
observed at kilometer wavelengths. Only the long wavelength elements of the 
herringbone structure are shown (from Cane et al., 1931). 

Figure 2. Histograms of the distribution of flare longitudes of the flares 
associated with the I? type II bursts. 

Figure 3. Count rate of the >13 Mev/nucleon detector on ISEE-3 (courtesy of 
T. T. von Rosenvinge). occurrences of IP type II bursts and SA events are 
indicated . 

Figxire 4. Histogram of the distribution of transit velocities of the shocks 
associated with the X? type II events. 

figure 5. Shock transit velocities shown as a flmction of the longitude of the 


associated flare 
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The long period variations of the first eight planets 
in the solar system are studied. First, the Lagrangian 
solution is calculated and then the 

long period terms with fourth order e^entricities and inclin- 
ations are introduced into the perturbation function. A sec- 
ond approximation was made taking into account the short 
period terms' contribution, namely the perturbations of 
first order with respect to the masses. Special attention 
was paid to the determination of the integration constants. 

The relative importance of the different contributions 
is shown. It is useless for example to introduce the long 
period terms of fifth order if no account has been taken of 
the short period terms. Meanwhile, the terms that have been 
neglected would not introduce large changes in the integration 
constants. Even so, the calculation should be repeated with 
higher order short period terms and fifth order long periods. 
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LOHC PERIODIC TERMS IN THE SOLAR SYS*r*E?! 

P. Bretagnon 

Bjreau des Longitudes, Observatoire 
17, rue E. Oeutsch de la Heurthe 
F-75014 Paris, France 

Received June 27, 1973 

SUMMARY [English language sunanary from the original text] 

We have studied the long period variations of the eight planets 
of the solar system (Pluto is excluded) . We first calculated the 
Lagrange solution. We then introduced the long period terms of 
fourth order in excentricities and inclinations in the disturbing 
function. In a second approximation we took into account the contrib- 
ution of the short period terms which provide the perturbations of 
the first order with respect to the masses. We have paid special 
attention to the problem of the determination of the integration 
cons tants. 

We began with the expansion of the disturbing function R [form- 
ula (l)J. We used the variables h = e sin a, K = cos p = sin ^ 

cos n, q = sin ^ cos fi and obtained expression (3) for the disturbing 
function and the equations of Lagrange (4). 

In the Lagrange method, one retains only the second order terms 

of the quantities h, k, p, q of the so called long period part of 

the disturbing function. The resolution of the system of differential 
equations thus obtained gives the solution of Lagrange (5) . The 
corresponding integration constants are given in Tables 2, 3, 4 and 5. 

We later introduced the long period terms of the disturbing 

function, of fourth order in the quantities h, k, p, q. These terms 

give rise to third order terms in the Eq. (6) for the variables h^, 
for example. We then substitute rumerically the Lagrange solution 

♦Numbers in the margin indicate pagination in the foreign text. 
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in these thiri order teras and hence obtain the form (8) of the 1 

equation for dh^/dt. 

In a second approximation, we also introduced the short period 
tenss of the disturbing function. The masses are substituted nuE?er- 
ically and the terns thus found are indentical in form to those 

I 

arising frcsn long period terms of fourth order of the disturbing i 

function and are directly added to the Eq. (8) . ■ 


To solve the syst^s of Eq. (8) and (9) , we used the Xrylov- 
Bogolioubov method, which consists in seeking a solution of the 
form Oil) with a modification of the frequencies given by (12) . 
Through (12) and derivation of (11) we obtain (13). In addition, the 
substitution of (11) into (8) and (9) gives (14), so that we get the 
two expressions (13) and (14) for dh^/dt and dk^/dt; their third 
order parts are given in (15) by identification. It is then possible 
to determine the quantities M , N . B. and C. introduced in 
(11) and (12). 


} 


The solutions are given by (16) and (17) and in Tables 8 to 13. 

i 

The comparison between Tables 3 and 8 shows that the integration i 

constants have been greatly modified, particularly for the planets 

Mercury, Venus, Earth and Mars. This is due to the importance of 

third order terms for these planets. Table 9 gives the modifications 

and of the frequencies as well as the new values of these 

frequencies: g. =g. +B.; s. =s. +C.. Tables 10 and 11 show the 

amplitude of the Lagrange solution calculated with the new constants; 

Tables 12 and 13 show the amplitudes M and N n of the argu- 

u,.,“ u,.,6 ^ 

ments of higher order. 


This work displays the relative importance of the different 
contributions: it is, for example, Uoeless to introduce the long 

period terms of fifth order if one has not taken into account the 
short period terms. We have included the major contributions; the 
neglected terms would not introduce large modifications of the 
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constants of integration. However, the calculation should be repeated 
including long period terms of fifth order and short period terms of 
higher order. 


Key words: planetary theory, secular ?>erturbations 


There have been several studies of long period terms in the solar 
system* Stockwell, Harzer (1895), Hill (1897), and more recently 
Brouwer and van Woerkcra (1950) and Anolik et al. (1969). 

Brouwer and van Woerkom calculated the Lagrange solution for the 
eight planets and in particular investigated the Jupiter-Saturn case. 
This was a continuation of the work of Hill, who had determined a 
mean perturbation function on the basis of Le Verrier's findings. 
Brouwer and van Woerkom used this perturbation function, which had 
been extended to sixth order excentricities and inclinations, for 
Jupiter-Saturn. It is difficult, however, to determine the accuracy 
of their result because Hill empirically established some of the co- 
efficients . 

Anolik et al. dealt with the eight planet case by introducing 
all the perturbation function's long period terms up to fourth order 
excentricities and inclinations. 

Our goal was to evaluate the significance of the various long 
period terms according to their origin. We too dealt with only the 
eight planet problem. Pluto was neglected for several reasons. First 
of all, the generally accepted mass of Pluto, which previously had 
been 1/360,000 the solar mass, is now 1/1,8000,000 with a large 
uncertainty: 

= iSOOOOOt 600000. 

»'r 

Moreover, Pluto's radius vector can be less than Neptune's, with the 
result that expansions in a, the ratio of senimajor axes, of the 
perturbation function, are no longer convergent. Finally, the intro- 
duction of Pluto's influence causes the appearance of very large 
resonances between Neptune and Pluto whose physical character is 




unclear. 
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Lastly, we calculated the eight planet Lagrangian solution and 
then introduced the perturbation function's fourth order terns as 
well as the contribution of the hort period terms of first order 
with respect to the inasses. In addition, we particulary concentrated 
on the problem of determining the integration constants because of 
the significance of the terms modifying the Lagrangian solution. 


The expansions of perturbation function R that we used are those 
constructed by Chapront at the Bureau des Longitudes. They take the 
fomr of analytical expansions in e and sin i/2, where e represents 
the excentricity and i the inclination of the orbital plane relative 
to the plane of origin. 


R = 


V X jsin ] X (sin -? ] cos 


( 1 ) 


with 


=Ji 'i +y,Q, . 


X being the planet's longitude, w the argument of the perihelion, and 

the argtunent of the node. The subscript I refers to the inside 

planet and E to the outside one. The perturbation function's long 

period portion is that oart for which X and X_ are absent, i.e. in 

I L 


which 


= j 2 = 0. The summation with respect to the small quantities 


e^, Cg, sin ij/2' ^E/2 done starting with zero order terms 

and then 2 , 3 , . . . 


The orbit's descriptive elements (the semimajor axis a, the 
excentricity e, the inclination i, the node argument fi, and the 
perihelion argument “ are those of Newcomb. These elements are 
expressed relative to the 1850.0 ecliptic averaged over short periods, 
which will serve as our point of departure (t = 0 for 1850.0) in 
determining the integration constants of the sought after solutions. 
Also, we used more recent values for the masses of Venus, Earth, 

Mars, and Saturn than the ones Newcomb used. 


The mean motions n^^ are the average observed values. The semi- 
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major axes a, are related to the values of n. by the expression 
"l^l “ constant. Now, we need for the mean motions values from 
which the secular perturbations have been removed. We therefore 
calculated the secular perturbations on the basis of Chapront's 
and Simon's work concerning the construction of planetary theory with 
secular terms. 

In the end, we used for each planet the value n of the mean 
motion defined by: 

n = - 5n 

and the value a of the semimajor axis obtained by 

2 3 ^ ^ 

n a = constant. 

We have assembled the elements adopted for the eight planets in 
Table 1. 


Table 1 

Planetary Elements for 1850.0 



a, r/yr) 

a r/yr) 

nj (AU) 

n (AU) 

t* 

n 

i 

Q 

fit m 


538I0I6.5S93 

5381023.1732 

0.3S709867I3 

0.387098.3460 

0.20.360.3% 

75 07 I9:.37 

7 0007:00 

46 33 I2.*24 

60(K|OII) 


2t06 641.4171 

2 l0665l.*63l 

0.72.3.3322169 

0.7233298487 

0.00684458 

l29 27.34.5 

.3 2.3 35.26 

75 1947.41 

4O8.M.0 

Earth 1295977.4496 

1295975.6094 

1.000000021 

1.000000968 

0.01677126 

lOO 2l .36..30 

0 


3289(0 

Min 

689050.9354 

6S9059.;'<|7 

1.52.36914428 

1.52.36791.387 

0.09326685 

.3.3.3 17 52.37 

1 5102.42 

•■48 24 03.40 

3099000 


109256.6.3954 

10926.3.05033 

5.202S03945 

5.202600424 

O.IM825382 

1 1 54 26.72 

1 1841.81 

98 55 58.16 

104*355 


43996.20414 

43865.62112 

9.53X84.3653 

9.554827.367 

0.056068)75 

90 06 39.53 

2 29 39.26 

112 20 51.38 

349S 

'Uanu5 

I5426.092S 

1 f *'5 1 

I9.1822SI85 

!9.2i7io6i3 

0.ia69055 

:68 15 46.9 

0 4f) 20.54 

7.3 14IJ8.0 

22869 

tiquune 

7864.698 

784.3.328 

.30.05‘'342 

.30.111791 

0.0085082 

43 19 43.7 

I 4701.81 

1.30 08 00.2 

I“3l4 


[Commas in tabulated material are equivalent to decimal points.] 

We chose the following variables to analyze our problem: 

h-csinn. p=sin{sinO, (2) 

k = ecosn. g=sin]cosf2. 

This choice was made in order to avoid the appearance of quantities 
expressed in e and i in the denominators of the Lagrangian eouations. 

Such quantities could cancel each other out. In addition, this is 
necessary for the resolving process because in this wav the solutions /143 
are expressed formally through the use of these variables and, in the 
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algorithm of solution's construction, the second members always re- 
tain the scune polynomial form. 


The change in variaoles defined by (2) yields in the perturbation 
function in form (1) an expression of the form: 

cosOV/ + ( 3 ) 

where the summation is extended to such exponential values that 

+ S 2 + tj^ + t 2 + u^ + U 2 , where is the order at 
which it is desired to limit the calculations. 


For the variables defined in (2) the Lagrange equations are 
written: 


Jb 

Jf~ 

7i 

dt 

J<L 

dt 


0 

-<?-»* - 

cR 

/>(! - t-'l 

1 2 

cR 

kp 

cR 

4. 

kq 



nar 

ck r.ii~ 

[Tr,i“ 

-T 

^ 

{/. 



Wn-e-)'"' 

N 


(I-f-)' 

- <'R 

1 

1 

•4 

cR 

IfP 

cR 

l_[j 

dR 


na~ 

cli -^(1 


7 71 

2 n«i'< 1 - f*)' 

* 

2»?<7‘( t — e'l* 

2 cq 


I 

^R 

P 

- 

vR 


cR 

4. 

ph 

cR- 

Ti 

fti*(l - e 

77 ■ 

1 - 

>“T'V 

c/. 

2 »ui’( 1 - ^ 

ch 

■2»u-7i -e')' 

~c'k 


1 

(R 


<1 

cR 

qk 

cR 

+ 

cR 


4na'{ \ - 

2 

cp 

2 ncj-( 

1 -f-l' 

‘ * f/. 

2 mj-(l -e*l‘ 

Ph 

2 na^(l - 

2 dk 


1 da 2 

cR 

a dt iia^ 

C/. 

d). 

2 cR 

— n — 
dt 

na ca 


_ [ \-e 

»u-[l + (l 


')'■ i, <*R 1 / rR . cR\ 

— ^ ~'T I — ’.7 ^’ 7 f^ P ~ — ~ — I 

ih < k j 3»iJ*(l -e‘) '■ \ cp cq j 


(4) 


2 2 2 
where e = h + k . 


LAGRANGIAN METHOD 


For a planet of subscript u perturbed by the seven other planets 
of subscript v, the perturbation function is written: 


r <« I > M 


The first summation is extended to the planets inside the one under 
consideration, the second to the planets outside. We use the follow- 
ing notation: 



a. /A (A 


distance of the two planets) 




and 
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Limited to the second order, has the following expression; 

K = + ‘'I + -t id + PI + «r) 

+ BjkX + W + -f /»«P,> 

where C , A , B are functions of a = a /a , which is constant 
uv uv uv uv u V 

here. We thus obtain, with the notation; 




if 

I + w. 


r , njn, . _ 

[u. r]= if i<ti. 

1+MI, 

-ir = + I 

(it , 

= - I [u. v]i2 A + 

-^ = - X [u.v](2A^,q,-2A„q,), 

-^ = + I 

at ^ 


This system is written in matrix form as; 

m dK c u 

—~=ExK, — 7 — = -fxW, 

dt dt 


-dT=‘^^ 


dt 


= -txP, 


where H is the column vector with components hy 


i • • • f 




the column vector (k„ , k„, ...» k„ ) , P the column vector (o„ , p„. 

Me V N ■ Me V 

...» P„) » and Q the column vector (q„ , q„, ..., q„) . The sub- 
N ne V N 

scripts Me, V, ..., N represent Mercury, Venus, ..., Neptune, 
respectively. E and I are the matrices of the linear systems in 
excentricities and inclinations respectively. 


The conventional resolution of the two Lagrangian systems gives 
the eigenvalues; 

g^, i = 1, 2, ...» 8 for the excentricities; 
s^, i = 1, 2, ..., 8 for the inclinations. 
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One of the eigenvectors in the system of inclinations in zero. 
We will assume s^ = 0. 


We determine the eigenvalues X. . associated with g., and u.. 

1 ^ 3 ^ ^ 

associated with s^, which gives the Lagrangian solution: 

t 

hi= Y. 

j- 1 
8 

Z /.ijMjCos{gjt + p^) 
j‘ I 

8 

= V 
y- I 

8 

Z C<'S(.s/ + . 


P.= 2, Hij>ljS\n{Sjt + Sj] 


(5) 


Lastly, we calculate the 32 constants of integration , 8^ , 
Nj, from the values of h, k, p, q ?t t = 0. • 


We have assembled the eigenvalues g and s in Table 2, and 

in Table 3 we show the 32 constants of integration M, B, N, 6. 

Lastly, Table 4 gives the amplitudes of the Lagrangian solution roul- 
8 8 

tiplied by 10 : A..M. x 10 for the excentricities , and similarly 

^8 . . 
in Table 5, y N. x 10 for the inclinations. Frequencies g and s 

i j J 

are expressed in seconds per year. X.., u . . , M., N. are dimensionless 

ih 13 3 3 

numbers. 


Table 2 

Frequencies in "/yr 
(Lagrangian Solution) 



.. 


} 



1 

I 


+ 5.4fi|369 
+ 7.346581 
+ 17.531295 
+ 18.004584 
.1711401 
_ 22,286552 
+ 2.701787 
+ 0.6.43116 




5 199958 
6.5*1.487 
ly*462()5 
r.6.46m 
0 

25 ‘-!l'6 
:.904526 
o.^-:520 


Table 3 

Constants of Integration 
(Lagrangian Solution) 


j 

M 

P 

.V 

A 

1 

0.1814)040 

87 )) )i:37 

0.06274851 

18 15 28:76 

2 

0.0)9097)2 

)92 40 29.79 

0.00506380 

4)6 20 04.8) 

J 

0.0)056860 

3.42 56 5).67 

0.01222166 

2.54 17 05.60 

4 

0.07.40040J 

.4)6 06 34.8) 

0.025)99)8 

295 40 38.75 

5 

0.044)9426 

27 48 .44.76 

0.0)38.4974 

106 08 44.50 

6 

0.04S47744 

) .77 5) 04.46 

<l.t)0*.Nh789 

l26 19 41.4.4 

7 

0.0.4)40786 

106 1925.59 

O.OOSS02K6 

3)3 42.55.4) 

8 

0.0092.4780 

66 09 18.45 

0.005 88386 

201 0052.94 


(Commas in tabulated material are equivalent to decimal points.] 
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Table 4 

g 

X 10 . Amplitudes of Lagrangian Solution 


J 

i 

1 

2 

> 

4 

5 

6 

7 

8 ’ 

% 

Mercury 

IS14ICUO 

-23ISI06 

I55K97 

- P0495 

2414071 

1I3S4 

62612 

727- 

Venus 

63IS73 

1900712 

-1275337 

1497153 

1624499 

- 55447 

61654 

1105 

Earth 

4()5f<')5 

|4904t.>‘' 

1056860 

- 1492703 

1624452 

24'593 

65295 

1284 

Mar.\ 

6ft:s3 

264664 

3017351 

7300403 

1870379 

I6T034 

86626 

2053 

Jupiter 

703 

- 1055 

95 

- 55 

4319426 

- ; 563I.UO 

218379 

5968 

Sulurn 

627 

1088 

750 

840 

3404356 

4'37743 

199229 

6'36'' 

Uranus 

271 

265 

44 

46 

- 4384597 

- 181906 

3140786 

141054 

Neptune 

4 

10 

3 

3 

160418 

- 13561 

- 338902 

923740 


U . . N . X 10®. 
ID D 


Table 5 

Amplitudes of Lagrangian Solution 


X. J 
i N. 

■ 

2 

3 

4 

5 

6 

7 

8 ' 

Mercury 

6274851 

- 1781583 

204668 

58171 

1383974 

13920 

- 166549 

- 72367 

Venus 

591896 

506380 

-1341594 

- 343391 

1383974 

6022 

- 95883 

- 66215 

Earth 

426404 

408232 

1222166 

226117 

1383974 

140699 

- 86614 

- 64905 

Mars 

90534 

90894 

- 1794150 

2519918 

1383974 

452917 

- 62850 

- 61481 

Jupiter 

- 1038 

655 

9 

88 

1383974 

-3t5S7g 

- 47877 

- 58499 

Saturn 

- 1328 

925 

241 

9l6 

1383974 

T66‘99 

- 39034 

- 563SI 

Uranus 

1112 

477 

20 

86 

1383974 

- 34-90 

880286 

M7U 

Neptune 

28 

27 

2 

9 

1383974 

- 3858 

- 103566 

588354 


INTRODUCTION OF HIGHER ORDER TERMS 

We are now going to introduce the perturbation function's long 
period terms of fourth order h, k, p, q, as well as the perturbation 
function's short order terms. 


Fourth Order Long Period Terms 


By differentiation, these terms yield third order terms, and the 
Lagrangian equation for variable h^, for example, then has the follow- 
ing form: 


th. 


Ji 


y [m. « fc. - Py,0K- Pu- r. . <?. >} 


( 6 ) 
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where is a homogeneous third degree polynomial. 

Into polynomial we substitute the Lagrangian solution (5), 
whose numerical values are given in Tables 2, Z, 4, and 5: 


■ s 

k.= X X /,^A/^cosVj. 

• » 

= X J sin Oj . 4* " X ' 

where we have made , . = g . t + 3 . and 6 s . t + 5 . . Therefore only 

3 3 3 j 3 3 

the numerical values of amplitudes > .M . and ,j .N. appear in this 

U3 j ^33 

calculation . 

Among the values of the i and j subscripts of arguments 
(i = 1, 2, 8) and (j = 1, 2, 8), such a substitul^ion 

makes combinations appear in which at most only three values of sub- 
scripts i and j are involved. For example, there will be combinations 
of the type ^^'^5 " "^6^ * 

The expression X therefore has the form; 

I • « 

X + + (7) 

••• fi.h '* 

> ji-Ji it 

where C • j ■ . is a numerical coefficient. 

u,ii. . .ig, . . jg 

The summation over integers i and j is such that; 

i « 

I l'«l + X l/«l = • or 3 

••• m., 

We will designate that; 


and hence eauation (7) takes on the form; 


‘X [«.rlP.r= X^.A.*^- 

?• V.» 


to^(V- f 'l. 


We also make; 


10 
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£= + I if 
£ = - 1 if 



■ 



1 


I y 




Equation 


(6) is then written: 



l»« v.( 


( 8 ) 


Substituting the Lagrangian solution into the equation in 


dk^/dt similarly yields; 

dk 


V,‘ =- I 


~ X ^ ^ 


Z. ‘ 


sinl^. ft) . 


(?) 


We also calculate: 


dPu 

«f 

d(J: 

dt 


‘ ’ •* *.» 




V.9 


Short Period Terms 


Substituting the Lagrangian solution into the short period part 
of the perturbation function yields only short period terms that are 
first order with respect to the masses. It is only with the second 
mass order that we come across long period terms again. 

This time we have to consider for each planet the complete 
system of Lagrange equations (4) , which we will write for a planet of 
subscript u in the form; 


/146 
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(10) 


We determine the short period effects argument by argument. For 

a short period argument i>.^ + j being given integers, the 

functions F have the form; 

f, cos (/A, +7/^) + f , sin (IV. » 

where F and F are polynomials in h , k , q , h . k , p , q , whose 
C S JTJ uu uvv v v 

coefficients are functions of = a^/a^, n^ and n^. (In the special 

case in which one of the two integers i, j is zero, i.e., in the case 
in which the short period argument takes on the form i^^» the 
functions F and F depend on h , k , p , q , n for v = 1, 2, ..., 8 

and on the seven quantities = a^/a^ for u ^ v.) 

We therefore substitute the Lagrangian solution into equations 
(10) , which after integration yield a short period increase in the 
Lagrangian solution. Then by doing a Taylor expansion of the second 
members of equations (10) , we obtain second order terms vrith respect 
to mass after substituting the first order that we have just found. 

We will retain only the second terms' long period parts. 

Since the masses are always substituted for numerically, the 
terms thus found in the second members of the Lagrangian equations 
have the same form as those coming directly from the perturbation 
function's fourth order long period terms. 

In contrast to the case of the perturbation function's long 
periods, for which we kept the fourth period terms, the criterion 
for choosing short period terms is numerical. What we did was to 
retain the beginning of the h, k, p, q expansion of all the arguments 


1 •> 


V 
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causing changes in the Lagrangian solution frequencies of more than 
lo”^ "/yr. 


RESOLUTION OF THE SYSTE.MS OF DIFFERENTIAL EQUATIONS 


We saw that the contribution of the short period terms took on 
the same form as the terms coming directly from the perturbation 
function. We therefore have to resolve a system of differential 
equations having the form: 


Jt 

dt 


X [«.p]{2/4«A + B..Ic,}+ X-\., .cos(v.0). 

(8) 

r*« K-* 

(9) 


r*n i.# 


as well as a similar system for variables p and q . 


For that, we are going to use the Krylov-Bogolyubov method. This 
method consists of finding a solution of the form: 


with 


8 

i»i V.8 



I 

I 

■ 

I 

i« I 


V.# 

/i,yiV^cos<?j+ X 
v.» 




= s, + C,. 


( 11 ) 
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( 12 ) 


By differentiating system (11) and taking account of (12) , 
we obtain: 


13 


where 
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djt 

ill 

dk 


Jf = I '■uJ^ltCOSVjX{g. + Bj)+ 

I , 

H 

Jl ~~ Y. Bj)— 


s.* 


ig. s) - i^g^ + Ijg, + • •• 4 . +y, 5 , + ••■ + . 


Furthermore, equations (8) and (9) yield, by plugging in (11) 

»//!„ V- r f f ^ 1 

-7r=L X /,;A/;COSVV+ Y, K.w.*COS(V. 0 ) 

I l,= l V# J 


d( 


I 11 

+ B,, I y ;.,^A/.cosv>+ X *v/;.,,,cos(g-.r») [ + X v 

V.* JJ V.(t ■ 

“ = - X X '^-.r'^sinv’y+ X eSin(v.O) 


+ B. 


8 

J=l 


V.» 


I v.» 


Hence, we have two expressions, (13) and (14), for dh /dk and 
dku/dt. We make equal their parts that are of third order with 
respect to variables h, k, p, and q: 


X ^ gcos{^\0) 

y=i v.» 

= X [“•*■] X + B,, ^ Af; ,, gCOS(l,’.0)l 4 X V «cos(v\fl» 

**“ 1 I* ».» J V.« 

» 

X Bj/.^jMj$iny.~ X (ff- 


8 

V 

J~l 


K.t 


- X i- ■*«. X •'^y.v.»S'"(V’.f^) + fiy, X ^ ,sin(v- B). 

' *■' 1 «. • V.8 I ».e 


The method now consists of establishing argument by argument 
identities within each order. Two cases arise: 

1) The case in which the argument ('., 9 ) is equal to ij; ^ , i.e. 

8 8 

X I'-K X ly-i = >• 

» * I »i » 1 


( 13 ) 


(14) 


(15) 


we 


have: 
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I 


Here, we once again come across the Lagrangian solution argu- 
ments, and establishing identities between coefficients makes it 
possible to determine the new frequency values. The solutions then 
are expressed in Fourier series of these new arguments. Establishing 
the identities yields: 


and 




= X [«• •' ](2-4., A/. + By, • 


For a given '<1/^, subtraction of these two equations furnishes: 

\r»» ' 


The fact that 


and [u,v]B are not zero means 
‘ uv 


that M = M' , whatever u and j are. 

U,4^j u,ii-j 

We can then write: 

4r*s * r»« 

In the first member of this expression, we once again come across 

matrix E of the Lagrangian system. Subtracting the eigenvalue g^ 

from the principal diagonal means that the M values (u = 1, 2, 

U ^ V' j 

..., 8) will not be independent. We then let: 

M. , =0 


This is an arbitrary step in the Krylov-Bogolyubov method. It 

reduces to changing variables over the integration constants M ^ , 

N-. The choice of M. =0 does not specify the solution but 
J :) ' ■ j 

imposes the choice of a certain type of expansion for the coefficients 
of arguments Vj* 


Having made this choice, we then have for each j 

system of eight equations in eight unkowns; B. and K 

n u , 4^ j 


fixed, a 
(u j) . 
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2 ) 


Case in which the argument is random, 

i i'j+ £ i/j=3. 


Establishing 

yields: 


argument by 

r»» 


• 1 mm I 

argument identities in 

+ Af, ^ — itf. 'lA/, ^ 0 


i.e. such that: 


equations (15) 


and 




By subtraction, we obtain M = whatever u and arqu- 

ment , G) is. 

We can then write: 




This time there is no arbitrary step and the resolution of the 
system of eight equations in eight unknowns gives for each argument 
(i;j,e) the eight values (u = 1, 2, 8). 


We therefore have expansions of and k^: 

Z /.;Af,sinv’y+ X Af,^ ^siniv.<'» 

j’t v.t 

*.= Z ^■^jMjCOsij'j+ VcA/„^ gC0S(v.<>) 

V.# 

And similarly we find: 

J‘ I y f 

v.» 


(16) 


(17) 
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RESULTS AND DETERMINATION OF THE CONSTANTS OF INTEGRATION 

In Tables 6 and 7 we give the integration constants and frequen- 
cies of the solutions obtained solely from the second order h, k, 
p, q long period terms. Comparison of tables 3 and 6 show how great 
the contribution of the perturbation function's fourth order long 
period terms is, especially for Mercury, Venus, Earth, and Mars. 

Table 7 contains the frequency modifications anc C^, as well as 

the frequencies' new values: g. =g. +B., s. =s. +C.. 

1 1 1 1 1 1 

Table 6 

Constants of Integration (According to the Solution Based 
on the Perturbation Function's Second and Fourth Order Long Periods) 


i 

M, 

A 

.V, 


1 

0.IS454.H67 

S3 .32 24-: 16 

0.05887664 

11 5l 44.62 

2 

0.0l8t>4:7({ 

191 26 15.44 

0.00323843 

30; 59 26.61 

3 

O.OI204IOI 

318 18 13.0.8 

0.0096’327 

33 14.44 

4 

0.0631 1073 

307 01 46.20 

0.03227762 

2*5 46 53.61 

5 

0.042974SS 

27 1743.13 

0,01.384057 

If* O' 34.95 

6 

0.048427S2 

127 29 49.50 

0.00786457 

125 49 11.16 

7 

0.03210686 

lOO 44 43.97 

0.00880119 

316 00 16.25 

8 

0,00932559 

64 54 28.20 

0.00592OS9 

201 11 07.71 


[Commas in tabulated material are equivalent to decimal points.] 


Table 7 

Modifications and New Frequencies in "/yr 
(According to the Solution Based on the Perturbation Function's 
Second and Fourth Order Long Periods) 


1 

0. 


c. 

}, 

I 

-0.258373 

* 5,202996 

-0.44.39^5 

- 5.643943 


-0,00*3'2l 

* 7..345860 

-0.22O51O 

- 6."91897 

3 

-0.1 300.32 

* 17.201263 

-0,l5rC>4 

- 18.898909 

4 

-0.16916S 

* 17.835416 

-0.23l.362 

- l'.>674'3 

5 

+ 0.018087 

* 3.7N4S8 

0 

0 

6 

+ 0.322115 

* 22.608667 

- 0.6O69O7 

- :6.348083 

*9 

+ 0.078.361 

+ 2.780148 

-0.08.3205 

r • 

1 

8 

*0.009180 

* 0.642296 

-0.009136 

- 0.686656 


[Commas in tabulated material are equivalent to decimal points.] 
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Table 8 

Constants of Integration (Complete Solution) 


1 


A 

s, 

1*. 

1 

0.17791613 

87 03 02;09 

0.059627*2 

|2 08 |7.|5 

2 

0.021 CM749 

193 35 04.97 

0.003 l53.vS 

.3i'5 I3 17.19 

3 

0.00988829 

319 43 16.64 

0.01002295 

24-< 01 59.07 

4 

0.06115173 

307 48 39.56 

0.031302*9 

2‘* 5606.16 

5 

0.04341616 

28 .30 11.60 

0.01.38.3939 

1(>6 i39 11.5* 

6 

0.04S 14727 

127 42 54.52 

0.00785328 

1 25 3.8 34.23 

7 

0.03126134 

114 46 31.5.' 

0.0O8S0OS8 

3l6 i* 55.94 

8 

0.00899181 

72 05 25.03 

0.0058SK06 

20: r 15.59 


[Commas in tabulated material are equivalent to decimal points.] 

Table 9 


Modifications 

and New 

Frequencies 

in Vyr 

(Complete Solution) 

\ 

1 

B, 

9, 

C. 

1, 

1 

-0.262290 

- 5.1990*'; 

-0.4 109*9 

- 5.6 k 9.3* 

2 

-0.000490 

+ 7..346091 

-0.199640 

- ^*71027 

3 

-0.110749 

+ 17.220546 

-0.08.3094 

- 15.829299 

4 

-0.147321 

+ 17.85726.3 

-0.18265S 

- 1 *.818769 1 

5 

+ 0.495804 

+ 4.207205 

0 

0 ‘ ’ 

6 

+ 3.930206 

+ 26.216758 

-0.525894 

- 26.267070 j 

7 

+ 0.363394 

+ 3.065181 

-0.095511 

- 2.999837 

8 

+ 0.034747 

+ 0.667863 

-0.0l.39li 

- 0.691431 ‘ 


[Commas in tabulated material are equivalent to decimal points.] 

Tables 8 and 9 give the constants of integration and the frequen- 
cies for the complete solutions, i.e. the solutions that take the 
short periods into consideration. By comparing tables 6 and 8, we 
can see that the integration constants are once more greatly altered. /151 
Comparison of Tables 7 and 9 show that while the short period terms 
hardly change the frequencies related to the inside planets, the g^ 
and g^ frequencies on the contrary are changed to a much greater 
extent by the short period terms than by the perturbation function's 
fourth order long period terms. 

The modification of the constants of integration originates in 
the magnitude of the nonlinear terms found in particular in the 


18 


expressions of the elements related to the inside planets. Of course, 
we began by calculating these terms by numerical substitution of the 
Lagrangian solution in Tables 4 and 5. We then determined an analytic 
form of the expressions found so as to calculate the new integration 
constants. With the help of this analytic form and by making a first 
order Taylor expansion about the first values of the integration con- 
stants, we obtained integration constants of sufficient accuracy after 
several iterations. 


g Table 10 

A .M. X 10°. Lagrangian Solution Amplitudes 
il D 



1 


3 

4 

5 

6 

7 

8 

MCTCUIV 

)779l6t3 

-2554951 

145862 

- 142815 

2426473 

11310 

62320 

708 

\rtus 

U 19702 

2 104749 

- 1 193243 

1254088 

1632845 

- 55183 

61366 

1076 

Earth 

39S077 

1642622 

988829 

-1250360 

1632798 

246415 

64990 

1250 

Man 

65007 

291694 

2823123 

6115173 

1879988 

1609341 

86222 

1998 

iupiler 

689 

1163 

89 

46 

4341616 

-1555604 

217360 

5809 

Saturn. 

615 

- t200 

702 

703 

3421845 

4814727 

198300 

6557 

I'ranut 

266 

293 

41 

38 

-4407122 

- 181041 

3126134 

137298 

Neptune 

3 

11 

3 

2 

161243 

- 13497 

• 337321 

899181 



Nj X 10®. 

Table 11 







Lagrangian Solution Amplitudes 




1 


3 

4 

5 

6 

7 

8 

Mercury 

5962772 

- 1109444 

167847 

■ 72261 

■ 1383939 

■ 13894 

166511 

- 72418 

Venus 

562458 

31553S 

-1100237 

- 426565 

1383939 

6010 

95861 

- 66262 

Earth 

405197 

254218 

1002295 

280886 

1383939 

140436 

86594 

- 64951 

Mars 

86031 

56602 

- 1471377 

3130279 

1383939 

482015 

62835 

- 61532 

Jupiter 

986 

408 

7 

109 

1383939 

-315287 

47866 

- 58541 

Saturn 

- 1262 

5‘’6 

197 

- 1138 

1383939 

785328 

39025 

- 56429 

If ran us 

1057 

297 

16 

107 

1383939 

- 34725 

880088 

54753 

Neptune 

27 

l6 

1 

12 

1383939 

- 3851 

103543 

588806 

Lastly, we 

give the 

totality 

of our 

solution 

in Tables 

8 to 

13. 


Hence, Table 8 contains the 32 integration constants. Table 9 gives 


the B, and frequency modifications as well as the frequencies' new 
values: = s^ + C^. The amplitudes of the Lagrangian 

solution corresponding to the new constants are given in Table 10 for 
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/153 


i 


■% 


/154 


i 

i 
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8 8 
.M. X 10 and in Table 11 for g . .N . x 10 . Finally, Tables 12 and 
Ij 3 3 

13 contain the amplitudes M q and N , of the higher order argu- 

Vi # ^ Vi ^ v* ' 

ments When computing these terms, we retained only the argu- 

-4 

ments whose amplitudes are higher than 10 for the planets Mercury, 
Venus, Earth, and Mars, and 10 ^ for Jupiter, Saturn, Uranus, and 
Neptune. The zeros found in Tables 12 and 13 are amplitudes less 
than the retained percisions. 


Table 12 


Vcunirni 

Mercury 

Venus 

Earth 

Mars 

Argument 

Mcrcur y 

Venus 

Earth 

Mars 

Vi 

iSl 

0 

0 

0 

V* - 0J - 04 

0 

0 

0 

- 177 


0 

0 

0 

- 199 

V’4 - 0J + 06 

0 

0 

0 

- 240 

+04 

0 

0 

0 

453 

V4 — Oj +04 

0 

- 536 

276 

6177 

V*-04 

0 

0 

0 

287 

V'4-204 

0 

0 

0 

176 

V» 

0 

981 

- 1064 

- 13379 

9>4 - 0« + 0» 

0 

0 

0 

370 

V» + 0, -#4 

0 

0 

0 

762 

V* 

- 120 

456 

- 503 

0 

Vs-V*-v, 

0 

0 

0 

134 

V'4 4 04 - 0» 

0 

0 

0 

- 354 

Vj-20, 

254 

0 

0 

0 

V’4 + 0J - 04 

- 519 

4126 

-3367 

- 14852 

Vs-0, -0, 

- H4 

0 

0 

0 

V'4+0J -0t 

0 

0 

0 

2ll 

V3'0, +0, 

102 

0 

0 

0 

V4 + V5-V'': 

0 

0 

0 

118 

r,-0, +0J 

- 275 

0 

0 

0 

-Vt - Vt. 

0 

0 

0 

- 286 

+04 

150 

0 

0 

0 

Vj-2v'« 

0 

465 

- 378 

- '26 

V, 

14547 

2918 

;;so 

1092 

V'j •- V ’4 - v« 

0 

0 

0 

- 133 

rs+0, -0, 

- 2172 

0 

0 

0 

Vi - V'4 ■*- v» 

0 

0 

0 

184 

fs + fl, -0, 

0 

- 216 

1 82 

303 

V'l -01 -1-04 

396 

0 

0 

0 

Vs + Vt-V? 

0 

0 

0 

105 

Vi 01 0j 

312 

0 

0 

0 

•^5" V'7 

- 118 

0 

0 

0 

v’i “ 0j -*■ 0,1 

- 237 

0 

0 

0 

V'i~V’j + 

0 

0 

0 

- 180 

V’ J “ ^3 ^ ^^4 

ri 

1 

2457 

- 2354 

8864 

-<■4-0, +04 

1644 

0 

0 

0 

Vi 04 + 0» 

0 

0 

0 

148 

+0, 

- 398 

0 

0 

0 

V'l 

0 

0 

0 

- 323 
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Table 12 


Argument 

Mcrcur.y 

Vfnus 

Earth 

Mars 

Vj + «4 -fl* 

0 

0 

0 

- 162 

Vi — ^4 

115 

- 901 

857 

- 1422 

Vi ■*■ 

0 

0 

0 

102 

Vi + Vi “ V’l 

0 

0 

0 

137 

Vi + V 4 - V* 
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0 
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0 
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CONCLUSION 


In this study of the long period variations of the olanetary 
elements, we added to the Lagrangian solution the terms of third 
order excentricity and inclination arising from the long period 
portion of the perturbation function calculated for the planets as 
a whole. We also took into consideration the influence of short 
period terms of second order mass. We particularly concentrated on 
determining the integration cc>nstants that make the solutions agree 
with the mean elements when t = 1850.0 is used as time zero. 

The terms calculated with these constants are grouped together 
in Tables 8 to 13. Notice in these results the very strong coupling 
that exists, for a long period problem, in the planetary system. The 
magnitude of the terms arising from the short periods shov;s that 
there is no point to extending a theory to the fifth order on the 
basis of the perturbation function's long periods if the short periods 
are not taken into account. 

This work's essential task was therefore the comparison of the 
various effects according tf their origin so as to have an overall 
view of this problem and to be able to embark on the complete construc- 
tion of a long period theory. Our solution is in fact still incom- 
plete. Even so, we should take into account the direct terms of fifth 
order that must have an influence, especially for Mercury , Venus, 

Earth, and Mars. We have yet to calculate the influence of short 
periods of higher orders of excentricity and inclination, and maybe 
even part of the third order with respect to masses in the case of 
the resonant argument 2X , - 5- between Juniter and Saturn. Such an 
investigation would be very important. Hov/ever, since the largest 
contributions have already been considered, it would no longer present 
any great difficulties for the determination of integration constants. 
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